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ABSTRACT 

The normalized change in ultrasonic "natural" velocity as a function of 
stress and temperature was measured in a unidirectional laminate of 
T300/5208 graphite/epoxy composite using a pulsed phase locked loop 
ultrasonic interferometer. These measurements were used together with the 
linear (second order) elastic moduli to calculate some of the nonlinear 
(third order) moduli of this material. 

INTRODUCTION 

Nonlinear elastic moduli are important physical properties in 
conventional materials. They provide information about the interatomic 
bonding forces in crystalline solids. Nonlinear properties are also 
important in the nondestructive determination of applied and residual 
stress (strain) and several investigations have also established a possible 
relationship between nonlinear elastic properties and ultimate strength in 
aluminum and carbon steel [1,2] . Conventional mechanical testing does not 
have sufficient sensitivity to measure the nonlinearity in the mechanical 
behavior in most materials. However, the effects of nonlinear elasticity 
on elastic wave propagation such as harmonic generation, elastic wave 
interaction, and velocity dependence on applied stress or temperature can 
be used to characterize these properties. 

In this research, an effort was made to determine the nonlinear elastic 
moduli of a unidirectional T300/5208 graphite/epoxy composite. This was 
accomplished by first determining the linear elastic moduli by both strain 
gage measurements during quasi-static uniaxial loading and ultrasonic 
velocity measurements. The stress and temperature dependence of the 
ultrasonic "natural" velocity was measured using a pulsed phase locked loop 
interferometer. Equations were derived to relate the stress dependence 
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of the velocity, and the linear elastic moduli to the nonlinear elastic 
moduli. In addition, a least squares method was applied to calculate a 
best fit for the nonlinear moduli from the measured data. Although 
unidirectional composites are most often characterized as transversely 
isotropic materials, the measured linear and nonlinear elastic properties 
exhibited small deviations from this transversely isotropic behavior. 

Thus, the nonlinear moduli were calculated for models of transversely 
isotropic and orthotropic elastic symmetry. 

LINEAR (SECOND ORDER) ELASTIC PROPERTIES 

The material used in this study was cut from a thick (0.8 in.) 
unidirectional laminate of T300/5208 graphite/epoxy. The original laminate 
was ultrasonically C-scanned for gross voids and defects. None were found 
to be present. 

The linear elastic stiffness moduli were calculated from measurements of 
the density and ultrasonic wave speeds at 2.25 MHz. Although these 
measurements would have ideally been made on a single specimen to avoid any 
errors due to sample to sample variations, the large attenuation of shear 
mode waves necessitated the use of several specimen. The dimensions and 
densities of these samples are given in Table 1. The longitudinal 
measurements were made on sample 1.10A while the shear wave measurements 
were made on the other three. The equations used for these calculations 
were similar to those presented by Kriz and Stinchcomb [3] for a 
transversely isotropic composite material. However, in this work, the 
fiber axis was designated as the x 3 direction as indicated by Fig. 1. The 

ultrasonic measurements demonstrated deviations from transversely isotropic 
behavior. Therefore, the stiffness moduli were also calculated using 
orthotropy as the elastic symmetry model. The stiffness moduli for this 
material for both the transversely isotropic and orthotropic symmetry 
models are presented in Table 2. The large uncertainties in the c 13 and 

c 23 moduli are due to the large uncertainties in the velocity measurements 
of the off axis non-pure mode waves needed to calculate these moduli. 

These measurements were made on specimens cut with the direction of 
propagation 45 degrees between the x 3 and x 3 axes. 

The linear elastic compliance moduli were determined using strain gage 
measurements and uniaxial compressive quasi-static loading. The measured 
moduli are presented in Table 3. Again, there was a small deviation from 
transversely isotropic behavior. The stiffness moduli were inverted and 
compared with the compliance moduli and were found to be in good agreement. 

NONLINEAR (THIRD ORDER) ELASTIC PROPERTIES 

Using the linear elastic stiffness and compliance moduli, the nonlinear 
elastic moduli can be calculated from measurements of the stress dependence 
of the ultrasonic "natural" velocity. The "natural" velocity (W) was 
defined by Thurston and Brugger [4] as the velocity referred to the 
unstressed or natural state. It is given by 

L 0 

w = — (1) 

where L 0 is the specimen length in the natural or unstressed state and t is 
the time of flight of the ultrasonic wave. Since L 0 is a constant, the 
normalized change of "natural" velocity with respect to a change in 



temperature or stress is given by 


Thus, only the change in time of flight of the ultrasonic wave need be 
measured . 

In these measurements however, a pulsed phase locked loop (P2L2) 
ultrasonic interferometer was used to measure changes in "natural" velocity 
at frequencies near 2.25 MHz. This instrument, developed by Heyman [5], is 
shown schematically in Fig. 2. The heart of the P2L2 is a voltage 
controlled oscillator (VCO) which generates a continuous wave signal of a 
frequency that is controlled by a D.C. input signal. A portion of the 
signal from the VCO is gated out into a tone burst and is then used to 
excite the ultrasonic transducer. The resulting elastic wave is launched 
into the sample and is then detected either by the same transducer in the 
pulse echo mode or another transducer in the through transmission mode. 

The received signal is input into the P2L2 where it is phase compared with 
the signal from the VCO at a preselected phase point using a sample and 
hold. The sampled voltage from the phase detector is then used to drive 
the VCO to a condition of quadrature. The acoustic phase shift (0) given 
by 

0 = 2 Jtf t (3) 

where f is the frequency, is then maintained as a constant by the feedback 
loop and therefore 

A0 = 0 = 2jt(tAf + fAt) . (4) 

Combining equations (3) and (4) yields 

A0 Af At 

— = 0 = + (5) 

A ft 


and thus. 


At AW 

t W 


By monitoring the normalized change of frequency of the P2L2, the 
normalized change in "natural" velocity is determined. 

The measurements made in this work were of the Stress Acoustic Constants 
(SAC) and the Thermal Acoustic Constants (TAC) . The Stress Acoustic 
Constant, as defined by Heyman [2] and Cantrell [6] is given by 


Aw Af 

W f 

H = = (7) 

Act Act 

where o is the applied stress which is either hydrostatic or uniaxial. The 
SAC's were measured for twenty seven combinations of stress, direction of 
propagation and polarization of the ultrasonic wave. To differentiate 
between these, each SAC is given three subscripts. The first represents 
the direction of applied stress which was limited to being along one of the 
three coordinate axes. A 0 for this subscript implies that the loading was 
hydrostatic. The second subscript gives the direction of propagation and 



the third gives the polarization direction. These were also limited to 
being pure mode waves along the three coordinate axes. The direction of 
propagation was always perpendicular to the direction of loading for the 
uniaxial SAC measurements which were made under compression. The same 
specimen that were previously described were used for these measurements. 

During the uniaxial SAC measurements, the specimen were loaded to 
between 80 and 100 MPa. The loading rate was varied considerably (25 to 
500 MPa/Min.) to determine its effect on the measured velocity changes. 
Within this range there was no measured effect indicating that viscoelastic 
effects were negligible. Each uniaxial SAC value was determined from the 
average of ten measurements during which the transducer was rebonded to the 
sample several times. Although most of the measurements demonstrated 
linear relationships between normalized frequency shift and stress, several 
exhibited nonlinear (quadratic) behavior. These were fit to a quadratic 
fit with the SAC value taken as the linear term of this fit. The measured 
values for the uniaxial SAC's are shown in Table 4 with those which were 
quadratic designated by an *. Typical curves for two of the uniaxial 
measurements are shown in Figs. 3 and 4. A block diagram of the 
experimental apparatus is shown in Fig. 5. 

The hydrostatic SAC measurements were made in a similar fashion except 
that the pressure chamber was limited to a maximum of 250 psi. This was 
about a factor of fifty lower than the uniaxial maximum stress. This much 
lower pressure produced smaller velocity changes making these measurements 
much more uncertain. This limitation in pressure also made the measurement 
of ^033 impossible because the changes involved were too small to be 
detected with the P2L2 . The measurements also had to be made over long 
periods of time to avoid temperature variations caused by increasing the 
pressure of the Argon gas in the chamber. Because of the long length of 
time required to make each measurement only three measurements were 
averaged for each hydrostatic SAC value. The results are presented in 
Table 4 with a typical curve shown in Fig. 6. The block diagram of the 
measurement system is shown in Fig. 7. For these small pressure changes, 
the relationships were all linear. 

The effect of temperature on the "natural" velocity was also measured 
for small variations about room temperature. The apparatus used for these 
measurements was the same as that used in the hydrostatic measurements. 

The temperature was increased above room temperature several degrees C in 
0.3 C steps using the heaters and temperature controller in the pressure 
vessel system. The frequency shift was measured at each point after 
thermal equilibrium was attained. The Thermal Acoustic Constants, defined 
as 

Af 

H = — (8) 

T AT 

where T is temperature, are given in Table 5. The second and third 
subscripts give the direction of propagation and polarization respectively. 
They were measured for nine combinations of propagation direction and 
polarization. The curves all exhibited linear relationships. 

The Stress Acoustic data was used to determine some of the nonlinear 
elastic moduli for this material. Since the linear elastic moduli and the 
SAC data exhibited small deviations from transversely isotropic behavior, 
the moduli were calculated for both the transversely isotropic and 
orthotropic elastic symmetry models. For the case of transverse isotropy. 



there are nine independent third order moduli while for orthotropy there 
are twenty. However, because Hg 33 was unable to be determined, several of 
the moduli could not be calculated. The transversely isotropic c 133 and 


c 333 and the orthotropic c 133 , c 233 , c 333 , c 123 


and c 456 could not be 


determined . 

The theoretical basis for these calculations comes from the work of 
Thurston and Brugger [4] who presented generalized equation for the 
derivative of the "natural" velocity with respect to stress evaluated at 
zero stress ( (p 0 W 2 ) p= g) . These equations were written in terms of the 

unstressed density (p 0 ) , the adiabatic linear elastic stiffness moduli 
( c S i jki ) , the isothermal linear elastic compliance moduli (s T i j kl ) , and the 
third order elastic stiffness moduli (c ijklmn ) . For the case of hydrostatic 


pressure the equation was 


(pi ) „ = 1 + 2wF + G 

p=0 HC HC 


( 9 ) 


where 


w = (p A W 2 ) n = C S N N U U , 

1 0 p=0 prqs p q r s 


F = s U U , 

HC aars r s 


( 10 ) 

( 11 ) 


and 


G r ,„ = S C N N U U . 

HC aauv uvprqs p q r s 


( 12 ) 


Np and U i are the direction cosines of the directions of polarization and 
particle displacements referred to the "natural" state. The summation 
convention on repeated indices is assumed throughout this paper. For 
uniaxial stress, the equation becomes 

- <Po w 2 >p=o = 2wF uc + G uc ' 


( 13 ) 



where 


F 


uc 


s T MMUU , 

abrs a b r s 


( 14 ) 


G =s T c MMNNUU , 

UC abuv uvprqs a b p q r s 


( 15 ) 


and M i are the components of a unit vector in the direction of stress which 
is always perpendicular to the direction of propagation of the ultrasonic 
wave. The relationship between these equations and the SAC is given by 


H = 


<PQ W >p=0 

(2p n W 2 ) 


( 16 ) 


Thus, a system of equations can be derived for the models of transverse 
isotropy and orthotropy that relate the measured linear elastic moduli and 
SAC's to the unknown third order elastic coefficients. These equations, 
although too lengthy for presentation here, were presented by Prosser [7] . 

A least squares procedure as outlined by Hankey and Schuele [8], was 
used to reduce the data and determine a best fit for the values of the 
nonlinear coefficients. This was done for both the transversely isotropic 
and orthotropic elastic symmetry models. The results of this reduction 
with the estimated maximum and probable limits of error are shown in Table 
6. The values are all negative with the exception of c 123 . This is in 

agreement with the nonlinear behavior of conventional materials. The large 
uncertainties in some of the moduli arise from the large uncertainty of the 
hydrostatic SAC measurements. These large uncertainties mask any 
deviations from transversely isotropic behavior in the third order elastic 
coefficients within the uncertainty of the calculations. 


SUMMARY AND CONCLUSIONS 


This research provides initial measurements of the characterization of the 


nonlinear elastic properties of graphite/epoxy composites using ultrasound. 



These will provide the basis for studying the relationships between 
nonlinear properties and more important engineering properties such as 
ultimate strength, residual strength after impact and fatigue loading, and 
fiber-matrix interfacial strength. These measurements also serve as a 
basis for a study into the ultrasonic measurement of residual stress 
(strain) in graphite/epoxy composites which may arise because of the 
mismatch in coefficients of thermal expansion of fiber and matrix. Thus, 
these measurements of the nonlinear properties may aid in the development 
of needed nondestructive evaluation techniques for graphite/epoxy 
composites. However, from a practical standpoint, the measured SAC's 
values are probably more useful than the calculated nonlinear moduli. This 
is because they are more easily obtained with smaller experimental 
uncertainty and can be measured in an individual specimen instead of 
calculated from measurements made in several specimen. Also, the SAC's are 
the quantities that are more useful in experimental attempts to measure 
applied and residual stress and may be important in nondestructive 
evaluation of ultimate strength of materials. 
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Fig. 1. Illustration of axis designation 
in a unidirectional composite 
with respect to fiber direction. 



Fig. 2. Block diagram of the pulsed phase locked 
loop ultrasonic interferometer. 






Fig. 5. Block diagram of uniaxial stress acoustic 
constant measurement apparatus. 



Fig. 7. Block diagram of hydrostatic stress acoustic 
and thermal acoustic constant apparatus. 


Table 1 - Sample dimensions and densities 


Sample Number Dimensions (in.) 



Density 


x i 

X 2 

x 3 


(g/ 

cm 3 ) 

1 . 10A 

0.8001 

0.7999 

0 . 8002 

1.5404 

+ / 

- 0.0004 

1 . 10B 

0.7999 

0.5002 

0.7996 

1.5460 

+ / 

- 0.0006 

1 . 8A 

0.8002 

0 . 8000 

0.5002 

1.5384 

+ / 

- 0.0005 

1 . 7 A 

0.5001 

0 . 8001 

0.7999 

1.5411 

+ / 

- 0.0005 


Table 2 - 

Linear Elastic 

Stiffness 

Moduli 



Transverse 

Isotropy Value 

(GPa) 

Orthotropy 

Value (GPa) 


c u 14.26 c 11 14.295 

C 12 6.78 c 12 6.78 

c i3 3.0 — 8.9 c-^2 


3.3 


9.0 







Table 2 


Continued 


Transverse 

Isotropy Value (GPa) 


Orthotropy 

Value (GPa) 

c 33 

108 . 4 


c 33 

108 . 4 

C 44 

5 .27 


C 44 

5.28 




C 55 

5.27 




C 66 

3 .74 




c 22 

14.226 




C 23 

6.6 - 7.7 


Table 3 - Linear 

elastic compliance moduli 

Transverse 

Value (GPa) -1 

Orthotropy 


Value (GPa) -1 

i — 1 
i — 1 

C/3 

0.089 

s n 


0 . 0890 

S 1 2 

-0.063 

s 12 


-0 .0626 

s 13 

-0.0021 

s 13 


-0 . 00208 

s 33 

0.00935 

s 33 


0 . 00935 



s 22 


0 . 0891 



s 23 


-0 . 00215 

Table 4 - 

Stress Acoustic Constants 



SAC Value (GPa) -1 SAC 

Value (GPa) -1 

SAC 

Value (GPa) -1 


H 122 * 0.0490 + /- 0.0009 H 211 * 0.0427 + /- 0.001 H 322 0.00116 +/- 0.00007 

H 311 0.00123 + /- 0.00001 H 121 0.0887 + /- 0.0007 H 212 0.0741 +/- 0.0010 

H 312 0.00279 + /- 0.00010 H 321 0.00299 +/- 0.00010 H 123 *0.068 +/- 0.002 

H 2i3* 0.0572 +/- 0.0010 H 323 -0.00993 +/- 0.0004 H 313 -0.00919 +/- 0.0003 

h 131 0.165 +/- 0.003 H 232 0.149 +/- 0.002 H 132 0.111 +/- 0.002 


H 231 


H 133 


H 233 


0.109 +/- 0.001 


0.0538 +/- 0.0005 


0.0479 + /- 0.0003 



0.39 +/- 0.02 


0.37 +/- 0.02 


0.257 +/- 0.007 


H 011 

H 021 0.261 +/- 0.006 

Hq3i 0.28 +/ — 0.02 


H 022 

H 013 0.27 +/- 0.03 
H 032 0.25 +/- 0.02 


H 012 

H 023 0.27 + /- 0.02 
H 033 Not measured 



TAC 


Table 5 - Thermal Acoustic Constants 
Value (1CT 4 C _1 ) TAC Value (10 “ 4 C” 1 ) TAC Value (10 “ 4 C _1 ) 


H T11 

-6.36 

+ /■ 

- 0.02 

H T22 

-6.33 +/- 

0.03 

H T12 

-8.5 +/- 

0 . 1 

H T21 

-8.5 

+ /- 

0 .01 

H T13 

- 10.2 +/- 

0.02 

H T23 

- 10.1 +/- 

0.2 

H T31 

-7 . 1 

+ /- 

0.5 

H T32 

-7.6 +/- 

0.7 

H T33 

-0.72 +/- 

0 .09 




Table 

6 - 

Third Order 

Elastic 

Moduli 




Transverse 

Isotropy 



Orthotropy 


Modulus 

Value 

(GPa) 

Maximum 

Error 

Probable 

Error 

Modulus 

Value 

(GPa) 

Maximum 

Error 

Probable 

Error 

c m 

-214 

17 

4 

c m 

-196 

15 

3 

c i 12 

-89 

12 

3 

c i 12 

-94 

5 

1 

C 1 1 3 

-4 

110 

23 

c i 13 

-63 

63 

13 

c 123 

65 

109 

23 

c 122 

-91 

8 

2 

C 144 

-33.4 

3 

0.5 

c 222 

-186 

21 

4 

C 155 

-49.1 

4 

00 

0 

c 2 2 3 

-60 

93 

19 

C 344 

-47 

30 

6 

C 144 

-33. 

0 2.4 

0.5 


c 244 -47.8 3.3 0.7 


C 344 -46 27 6 


c 155 -50.1 3.9 0.8 


C 255 —33.5 2.8 0.6 



c 355 

-49 

32 

7 

C 1 6 6 

-33.9 

0 . 6 

0 . 1 

c 2 66 

-33.1 

0 . 6 

0 . 1 

c 366 

-28.5 

4 . 0 

00 

o 




